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We use different semi-empirical dispersion correction schemes to assess the role of long-range 
van der Waals interactions in the adsorption of the prototypical molecular switch azobenzene 
(CeHs— N2-C6H5) at the coinage metal surfaces Cu(lll), Ag(lll) and Au(lll). Compared to pre- 
ceding density-functional theory results employing a semi-local exchange and correlation functional 
we obtain partly sizable changes of the computed adsorption geometry and energetics. The discom- 
forting scatter in the results provided by the different schemes is largely attributed to the unknown 
form of the damping function in the semi-empirical correction expression. Using the congeneric 
problem of the adsorption of benzene as a vehicle to connection with experiment, we cautiously 
conclude that the account of dispersive interactions at the metal surfaces provided by the various 
schemes is in the right ballpark, with the more recent, general schemes likely to overbind. 
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I. INTRODUCTION 

Since the advent of semiconductor based microelec- 
tronics, the component density of integrated circuits has 
increased at a steady exponential rate. As component 
sizes approach the nano-scale, maintaining this pace will 
require a shift towards alternative materials more ver- 
satile than doped silicon. One such approach, already 
on the research agenda for several years, is to exploit 
the great variety and adaptability of organic molecules 
in component design. In the context of such a molecu- 
lar nanotechnology, molecules with properties bi-stably 
and reversibly modifiable by external stimuli, so-called 
molecular switches, are a research topic of paramount 
importance. The azobenzene molecule (CgHs— N2— C6H5) 
qualifies in this class by undergoing a reversible trans-cis 
isomerization. Because of its relative chemical simplic- 
ity, this molecule has come under intensive experimental 
scrutiny. Numerous potential applications have been pro- 
posed, e.g. as a light-driven actuator^, or data storage 
medium^. On the theoretical side, the precise mecha- 
nism behind the switching from the planar, C2h symmet- 
ric trans isomer to the torsioned-twisted C2 cis isomer, 
is a much debated topic. The consensus view explains 
the observed trans-cis (cis-trans) isomerization, follow- 
ing photo-excitation of the tt — > tt* (n — > tt*) resonance, 
in terms of conical intersections between the ground and 
low lying excited states on the isomerization reaction 
coordinate^iii^. 

Nevertheless, gas phase or solvated switches carry their 
own set of limitations within the framework of molecular 
microelectronics, wherefore in recent years, properties of 
switches adsorbed at solid surfacesS' has emerged as an 
important research field. Since the switching function is 
an innate property of the azobenzene electronic struc- 
ture, the choice of substrate and substrate coupling is 
non-trivial: While bonding strong enough to localize and 
order switches is desirable, significant hybridization of 
the mentioned frontier orbitals, or steric hindrance due 



to substrate registry is not. Close-packed coinage metal 
(Cu, Ag, Au) surfaces ostensibly offer a reasonable such 
balance. However, even at these substrates, azobenzene 
and -derivate switches exhibit a host of modifications to 
the switching function, in all combinations of reversible 
or irreversible switching, and switching by light or exclu- 
sively by a local electric field^i^i^. 

In this situation first-principles calculations stand to 
offer a unique perspective on the geometric and elec- 
tronic structure of these systems. With such motivation 
we have recently performed a density-functional theory 
(DFT) study of the metastable states of azobenzene ad- 
sorbed at coinage metal (111) surfacesi^. Within the 
employed generalized-gradient approximation (GGA) to 
the DFT exchange-correlation (xc) functional, the essen- 
tial findings of this study are summarized as follows: At 
all three substrates the most stable adsorption geome- 
try corresponds to the azo (-N=N-) bridge centered on 
and aligned with the bridge site of the (111) lattice, cf. 
Fig. below. In the flat trans isomer geometry the 
two phenyl (— CgHs) moieties lie parallel to the surface, 
while in the three-dimensional cis isomer geometry they 
stand upright, pointing away from the surface. In this 
decomposition into azo-bridge and phenyl-ring moieties, 
the overall bonding is characterized by a balance between 
three major effects: The binding energy gained by the 
formation of a covalent-type bond between azo-bridge 
and surface, the energetic penalty due to distortion of 
the gas-phase molecular geometry upon adsorption, and 
Pauli repulsion between the phenyl-rings and the sub- 
strate. In the flat trans adsorption geometry the flrst and 
third effect are in conflict in the sense that the formation 
of covalent-type bonds requires shorter surface distances, 
where the phenyl-rings already suffer from strong Pauli 
repulsion. Since this conflict does not arise in the three- 
dimensional cis geometry, the latter is relatively more 
stabilized upon adsorption: While the gas phase trans 
isomer is more stable than cis by some 0.6 eV, this value 
is lowered upon adsorption and at Cu(lll) - where the 
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azo-bridge to surface bond is strongest ~ DFT-GGA even 
predicts a reordering of the isomer stability, with the cis 
isomer '--^ 0.3 eV lower in energj>i^. 

Within this picture of a delicate balance between com- 
peting effects, it is clear that these results will be sen- 
sitive to the well-known inabilitjsii of semi-local GGA 
xc-functionals to account for a fourth surface bonding 
contribution, namely dispersive van der Waals (vdW) 
interactions. For the present system a rough estimate 
of the importance of this contribution may be obtained 
from the similarity of the phenyl moiety and a ben- 
zene (CgHg) molecule. Comparing the ^^0.6 eV benzene 
binding energy derived from temperature programmed 
desorption (TPD) experiments— to the essentially 
zero value obtained within DFT-GGA^- , vide infra, sug- 
gests that the computed DFT-GGA binding energy of the 
trans azobenzene isomer with its two flat-lying phenyl- 
rings might be underestimated by more than 1 eV. While 
this underscores that a real understanding of azoben- 
zene and related organic molecules at metal surfaces 
can only be obtained from a properly balanced descrip- 
tion of all four surface bonding contributions, realizing 
this in electronic structure calculations for such systems 
is still a largely unresolved challenge. In fact, an ac- 
curate account of dispersion interactions is one of the 
major issues in contemporary large-scale first-principles 
modeling, and recent years have seen a plethora of pro- 
posed solutions, ranging from high-level wave- function 
techniques within various embedding schemes^^iii, via 
xc-functionals explicitly constructed to include disper- 
sion interactionsiSii^j to the electron correlation resulting 
from applying the random-phase approximation (RPA) 
in the context of the adiabatic-connection dissipation- 
fluctuation theorem^Sii^. 

Unfortunately, due to the sheer adsorbate-dictated sys- 
tem size none of these approaches are presently computa- 
tionally tractable for the problem at hand. In this situa- 
tion, much more modest semi-empirical schemes that cor- 
rect at least for the long-range vdW interaction represent 
a viable alternative for a first assessment of how much 
the missing fourth bonding contribution skews the previ- 
ously obtained DFT-GGA results. However, while com- 
putationally highly efficient, these semi-empirical disper- 
sion correction approaches22i2^iSiiS^i2£ are also not un- 
problematic. On the practical side their semi-empirical 
derivation has given rise to a manifold of suggested 
schemes that all have the same conceptual structure, 
but differ in their material-specific parameters. On the 
fundamental side the mere validity of these schemes for 
bonding at metal surfaces is uncertain. Apart from the 
system-specific interest a second objective of our study 
is therefore to subject the obtained dispersion corrected 
results to critical scrutiny. For this we analyze the scat- 
ter when applying different prevalent semi-empiricial cor- 
rection schemes to the azobenzene problem and use the 
adsorption of benzene at the close-packed coinage metal 
surfaces as a vehicle to isolate the balance between Pauli 
repulsion and vdW interaction within this methodology. 



The overall structure of the paper is correspondingly 
as follows: The next theory section provides a detailed 
description primarily of the different employed semi- 
empirical correction schemes and how they are integrated 
into the DFT calculations. Thereafter we present the 
changes obtained through these dispersion corrections to 
the previous pure DFT-GGA results for the geometric, 
energetic and electronic properties of adsorbed azoben- 
zene at Cu(lll), Ag(lll) and Au(lll). This is followed 
by a critical discussion which includes the detailed results 
when applying the semi-empirical correction schemes to 
the simpler problem of benzene at these three surfaces. 

II. THEORY 

A. Semi-empirical dispersion correction schemes 

In the semi-empirical dispersion correction approach, 
the missing dispersion contribution to the inter-atomic 
interaction is approximated by a simple isotropic poten- 
tial. At long range, this potential is given by the leading 
C'e.ij ■ Rij^ term of the London series, where Cg.ij is a 
material-specific, so-called dispersion coefficient between 
any atom pair i and j at distance Rij . At short range, the 
long-range expression is matched to the DFT potential 
by multiplication with a damping function f^R^^Rij), 
which reduces the additional dispersion contribution to 
zero, subject to a cutoff defined by some suitably cal- 
culated combination R% of the vdW radii of the atom 
pair, as illustrated in Fig. [TJ The dispersion corrected 
xc-functional is then formed by simply adding the cor- 
rection potential to the ordinary DFT xc-functional. As 
Ca^ij coefficients are additive^l, the dispersion corrected 
total energy i^tot may therefore be written as 

N N 

Etot = Euft+sq fiSnR^j, Rij) Ce^jj R^f , (1) 

i=l j>i 

where -Bdft is the standard DFT total energy, and the 
sums go over all N atoms in the system. In existing 
schemes, heteronuclear i?? and Ce^ij coefficients are ap- 
proximated from semi-empirically determined homonu- 
clear parameters, and differences between DFT xc- 
functionals in the description of short to medium range 
dispersion interaction are taken into account by a suit- 
able modification of the correction potential through the 
parameters sq or Sr as further detailed below. 

While effortless to evaluate computationally, this gen- 
eral expression has obvious weaknesses, of which perhaps 
the most glaring is the (free-) atom-in-molecule approxi- 
mation: The substantial variation of properties between 
effective 'atoms' of the same species in the molecule, e.g. 
the influence of hybridization states on effective polariz- 
ability, is neglected. This may be particularly severe for 
the metal surface, where the lowering of the effective di- 
electric constant due to screening^^ should be reflected by 
reduced dispersion coefficients of atoms in deeper layers 
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FIG. 1: (Color online) Illustration of the dispersion correc- 
tion between a Ag-C atom pair as a function of their in- 
teratomic distance Ra^-c- Shown in the lower panel is the 
S6C'6,Ag-c-RAg-c term in the OBS (medium dash-dotted line) 
and G06 (medium solid line) scheme, as well as the finally re- 
sulting correction after multiplication with the damping func- 
tion /(J?Ag-Ci Ras-c) (separately shown in the upper panel), 
cf. Eq. (HI: OBS (thick dash-dotted line) and G06 (thick solid 
line). Note the largely different form of the dispersion correc- 
tion, which at not too different Ce.Ag-c coefficients in the two 
schemes, cf. Table IH is predominantly due to the very shallow 
damping function employed in the OBS scheme. At distances 
comparable to the sum of interatomic vdW radii, 2.30 A and 
3.09 A for OBS and G06 respectively (marked as thin vertical 
lines in the upper panel), this leads to energetic corrections 
that differ by an order of magnitude. The TS correction po- 
tential is qualitatively very similar to that of G06, differing 
in seCe and SrR!^ parameters only. 



in the substrate. The apphcabihty of such semi-empirical 
dispersion corrections to (transition) metal bonding sys- 
tems is therefore uncertain. In order to qualify this 
to some extent for the present azobenzene at coinage 
metal surface problem, we therefore study the variation 
of our results between three published dispersion correc- 
tion schemes that (apart from different semi-empirical 
material parameters) each have their own approach to 
this issue: 



The first scheme is due to Ortmann, Bechstedt and 
Schmidt (OBS)^-, and employs the London formula^® to 
calculate the dispersion coefficients from experimentally 
measured polarizabilities and ionization potentials. The 
corresponding dispersion coefficients, Ce^ij are compiled 



TABLE I: Calculated Ca^ij dispersion coefficients and cut- 
off radii R^j between elements i and j. As the dispersion 
coefficients in the TS scheme depend slightly, ~ ±5 % for 
adsorbed trans azobenzene, on the effective atomic volume 
we only show averaged values over all atomic constituents of 
each species. To allow direct comparison between the effec- 
tive parameters employed by the three schemes, we show the 
values sgCs.ij and SnRij, i.e. after multiplication with the 
xc- functional dependent scaling factors (see text). Note that 
the G06 scheme does not feature parameters for Au. 





H C N 


H C N 


OBS 
Au G06 

(TS) 


31.9 77.7 54.0 
19.7 62.0 45.9 


1.81 2.21 2.19 
3.48 3.80 3.71 


OBS 
Ag G06 
(TS) 


35.0 86.1 59.2 
14.5 51.1 42.8 
19.3 60.9 43.7 


1.90 2.30 2.28 
2.64 3.09 3.04 
3.45 3.78 3.69 


OBS 
Cu C06 
(TS) 


30.6 75.0 51.6 
9.5 33.8 28.3 
17.9 55.6 40.1 


1.75 2.15 2.13 
2.56 3.01 2.96 
3.45 3.77 3.68 



in Table [H together with those of the other two schemes. 
The notably shallow, long-ranged damping function in 
the OBS scheme illustrated in Fig. [1] was fitted such 
that together with i^oFT within the GGA of Perdew and 
Wang (PWOl)"^" the experimental c-lattice constant of 
graphite is reproduced. Demonstrating the effect of such 
a shallow damping function is the main reason why this 
scheme is considered here, even though it was not devel- 
oped for universal transferability and does therefore e.g. 
not provide fits to other xc-functionals in the original 
publications^'^ (i.e. sg = Sr = 1). 

The second correction scheme employed is the 2006 re- 
vision proposed by Grimme (G06)^. It features param- 
eters calculated from first-principles for most of the pe- 
riodic table (although not including Au), a compared to 
OBS steeper and shorter ranged Fermi damping function, 
cf. Fig. [U and a xc-functional dependence determined by 
fitting to thermochemical benchmark calculations. This 
scheme does aim at transferability, and has garnered pop- 
ularity in the field, with applications including adsorp- 
tion problems^i^^ akin to this study. However, for the 
present work, two GOB design choices may prove detri- 
mental: First, the problem of defining transferable tran- 
sition element 'atoms', is admittedly^ crudely solved by 
taking their parameters as simply the average of those of 
the preceding noble gas atom and the following group III 
element. Second, the xc-functional dependence is intro- 
duced by (effectively) scaling the dispersion coefficients 
with a factor sq 1 {Su = 1), cL Eq. ([J), which corrects 
most for the xc-functional where its influence is smallest, 
i.e. it shifts the Cg^ij potential also at long range, 
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where it can be expected to be most accurate. 

The third scheme was recently put forth by Tkatchenko 
and Scheffler (TS)^, and represents the state of the 
art: Expfoiting the relationship between polarizability 
and volume^^, this scheme accounts to some degree for 
the relative variation in dispersion coefficients of differ- 
ently bonded atoms. This is achieved by weighting values 
taken from the high-quality first-principles database of 
Chu and Dalgarno^i, with atomic volumes derived from 
Hirshfeld partitioning^ of the self-consistent electronic 
density. The TS scheme uses the same damping function 
as the GOG scheme, but following Jurecka et al^, in- 
stead scales the combined vdW radius by a xc-functional 
dependent factor Sr ^ 1 (sg = 1), thereby correcting 
where the xc-functional influence is strongest, and leav- 
ing asymptotics intact. 



B. Density-Functional Theory Calculations 

All DFT calculations were performed with the 
CASTEP— code using a plane-wave basis together with 
ultrasoft pseudopotentials'^^ as provided in the default li- 
brary, and the GGA-PBE functional* to treat electronic 
exchange and correlation. For benzene adsorption the 
local-density approximation (LDA) in the parameteriza- 
tion by Perdew and Zunger— was additionally consid- 
ered. The computational parameters were exactly those 
used and detailed already in the preceding studyiS, which 
is why we only briefly recapitulate those aspects of rel- 
evance for the dispersion correction schemes here. For 
azobenzene and benzene we modeled the (111) surfaces 
as four layer slabs in supercells with a lateral period- 
icity of (6 X 3) and (3 x 3) surface unit-cells with 18 
and 9 atoms per layer, respectively. Dispersion correc- 
tions were implemented in an external module^^ such, 
that a boundary-condition dimensionality < D < 3 can 
be chosen for arbitrarily defined subsets of the supercell 
atoms. Aiming to describe properties in the low coverage 
limit, this allows to account for dispersion corrections be- 
tween the adsorbate and an extended substrate, while si- 
multaneously suppressing spurious long-range dispersion 
interactions between the adsorbate and its periodic im- 
ages. Where applicable, parameters optimized for the 
employed GGA-PBE functional were used in the dis- 
persion correction schemes, i.e. sg = 0.75 (G06}^ and 
Sr = 0.94 (TS)?5 ^as used in Eq. ©. 

The module was interfaced to a locally modified ver- 
sion of the academic release of CASTEP 4.3 allowing to 
include the externally calculated dispersion corrections 
to total energies, forces and stress. Furthermore, it was 
extended by routines providing the atomic volumes based 
on a Hirshfeld analysis"^^ as required for the TS scheme. 
The Hirshfeld analysis is based on the (soft) pseudo 
charge densities for both the superposition of atomic and 
the self-consistent densities. Hirshfeld charges and vol- 
umes are integrated on the (standard) real space grid 
as provided by the plane wave basis set. We carefully 



checked the convergence of the volumes with respect to 
the size of the simulation cell by integrating over appro- 
priate supercells (using properly replicated charge densi- 
ties in the self-consistent case and applying some special 
treatment for hydrogen atoms to suppress the numerical 
noise on the integration grid). For several small organic 
molecules and the fee bulk phases of Au, Ag and Cu we 
have compared our implementation to reference values 
obtained with the all-electron code FHI-AIMS"*^ which 
employs numerical atom centered basis functions based 
on radial grids. For the transition metal (pseudo-) atoms 
the obtained Hirshfeld volumes differ rather substantially 
from the reference values. However, the ratios of the two 
volumes, which are the only required input for the TS 
scheme, are in very good agreement with the reference 
numbers in all cases and for all elements in this study. 
As a side remark, we note that this demonstrates for the 
first time the applicability (with very modest implemen- 
tation effort) of the TS scheme to DFT calculations done 
with a plane wave basis set. Further details on this im- 
plementation and the numerical tests will be published 
elsewhere''^ . 

As in our preceding studj^ the centrally targeted en- 
ergetic quantities are the adsorption energy i?ads and the 
relative isomer stability AE. The former quantity is de- 
fined as 

Eads — 2 [^azo@(lll) — £-(111) ^ -E'azo(gas)] , (2) 

while the latter quantity is defined as 

AE = ^ [Sazo@(iii)(cis) - £;azo(Q)(iii) (trans)] . (3) 

Here, -Eazo@(iii) is the total energy of the relaxed 
azobenzene-surface system, i^(iii) the total energy of the 
clean slab, -E'azo(gas) the total energy of the correspond- 
ing relaxed gas-phase isomer (all three computed at the 
same plane- wave cutoff), and the factor 1/2 accounts for 
the fact that adsorption is at both sides of the slab. The 
adsorption energy of cither cis or trans isomer at the 
surface is thus measured relative to its stability in the 
gas-phase, and a negative sign indicates that adsorption 
is exothermic. AE > 0, on the other hand, indicates 
that the trans isomer is more stable at the surface. Sys- 
tematic convergence tests indicate that these quantities 
are converged to within ±30meV at the chosen compu- 
tational settings^". Defined as differences of surface with 
and without adsorbate resp. adsorbed trans and cis to- 
tal energies, i?ads and AE are virtually unaffected by 
the intra-substrate, i.e. metal-metal dispersion correc- 
tions. Residual contributions arise from the only slightly 
different surface atom relaxations upon adsorption and 
in case of the TS scheme from different Hirshfeld de- 
compositions. Given the anticipated limitations of the 
semi-empirical approach in describing metal surfaces, we 
deliberately switch these small contributions off, i.e. the 
reported dispersion-corrected i?ads and AE only consider 
intra-molecular and molecule-substrate vdW contribu- 
tions. Here, the intra-molecular dispersion corrections to 



these quantities arising from the changed molecular ge- 
ometry in the adsorption structure are in fact minute and 
are below 15meV at all surfaces and for all schemes. For 
benzene, we additionally calculated the zero-point energy 
correction to the adsorption energies E^g^. The required 
normal modes of the benzene molecule have been ob- 
tained as F-point phonons. In our supercell setup these 
are computed via finite displacements from equilibrium 
positions to obtain the nuclear Hessian by numerical dif- 
ferentiation of the resulting forces. 

III. RESULTS 

A. Gas-Phase Azobenzene 

Since the azobenzene adsorption properties are con- 
sistently referenced to the gas-phase molecule, it forms 
a natural starting point for our analysis. As previ- 
ously reportedi^, GGA-PBE yields the gas-phase ge- 
ometric structure of both isomers in close agreement 
with B3LYP^^ hybrid- functional DFT, many-body per- 
turbation theory calculations and experiment. This re- 
sult remains unchanged here; none of the correction 
schemes applied modifies the pure GGA-PBE gas-phase 
geometries significantly. The gas-phase cis-trans rela- 
tive energetic stability is also only modestly reduced, 
from AE = 0.57 eV for pure GGA-PBE, to 0.44 eV 
(OBS), 0.47 eV (G06) and 0.49 eV (TS) in the three semi- 
empirical schemes. This result is readily explained by the 
shorter distance between the phenyl-rings in the three- 
dimensional cis isomer, which leads to a larger dispersion 
stabilization in the cis isomer and thereby to a lowering of 
AE. This lowering by ~ 0.1 eV due to the approximate 
account of dispersive interactions suggests that the near- 
perfect agreement between GGA-PBE, AE = 0.57 eV, 
and experiment, AE = 0.6 eV—, is due to a fortuitous 
cancellation between long- and short-range errors in the 
semi-local functional. This interpretation is corroborated 
by the AE = 0.68 eV^ obtained with B3LYP, which is 
intended to improve on the DFT-GGA at short range. 
As the latter overshoots by ~ 0.1 eV, one could suspect 
that dispersion corrected B3LYP would then restore the 
agreement to experiment. 
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FIG. 2: (Color online) Top panel: Top view of the azoben- 
zene trans isomer at the preferred DFT-GGA adsorption site. 
Bottom panel: Side view of the corresponding cis isomer ad- 
sorption geometry. Additionally shown in both panels are 
definitions of analyzed geometry parameters z, (Inn, i^, and 
a (see text). 

using all three dispersion correction schemes, show a clear 
energetic preference for the previously determined most 
stable adsorption site at GGA-PBE level of theory. In 
this geometry, the central azo-bridge is aligned in a 1:1 
N-metal coordination as shown in Fig. [51 The resulting 
adsorption geometry is suitably described by four charac- 
teristic parameters also shown in Fig. [21 namely the ver- 
tical distance from the top surface layer to the azo-bridge 
plane z, the NN bond length cJnn , the CNNC dihedral an- 
gle w, and the CNN bond angle a. At DFT-GGA level 
the latter angle a remained unchanged from its gas-phase 
value, 115° and 124° for trans and cis azobenzene respec- 
tively, for both phenyl-rings at all three substrates. All 
three semi-empirical correction schemes similarly leave 
this angle unaffected, which is why Table [Til concentrates 
on the obtained results for the three remaining structural 
parameters. 



B. Adsorbate geometric structure 

Turning to azobenzene adsorbed at the coinage metal 
surfaces, we first note that due to the long-range nature 
of the dispersive interactions the energetic correction pro- 
vided by the semi-empirical schemes only shows a small 
corrugation with respect to lateral movements of the ad- 
sorbed molecule across the surface, leaving the optimal 
adsorption geometry largely determined by the DFT en- 
ergetics. Correspondingly, full geometry optimizations 
of all high-symmetry adsorption geometries as defined in 
Fig. 1 of our previous publication's, at all surfaces and 



At DFT-GGA level gas-phase molecular geometries at 
Au(lll) are largely unmodified by adsorption, a fact 
which dispersion corrections do not change: OBS does 
not affect the trans isomer at all, and the geometry cor- 
rection in the other cases is essentially limited to a small 
rigid shift of the adsorbate towards the substrate. For 
the trans isomer this means that the dihedral angle uj 
remains at its gas-phase value of zero, i.e. the molecule 
stays planar. At Ag(lll), the picture is somewhat more 
varied: Here, the stronger azo-bridge— surface interaction 
activates and elongates the NN bond in the cis isomer, a 



TABLE II: Azobenzene structural parameters as defined in 
Fig. [5] and as obtained using GGA-PBE, as well as the three 
semi-empirical correction schemes due to Ortmann, Bechst- 
edt and Schmidt (OBS), Grimme (G06) and Tkatchenko and 
Schefiler (TS). None of these schemes affects the gas-phase 
geometric parameters, which is why only the PBE values are 
quoted here. Note that the G06 scheme does not feature pa- 
rameters for An. 





Trans 

Z diMN 

(A) (A) (deg) 


Cis 

Z ciNN 

(A) (A) (deg) 


Gas-phase PBE 


- 1.30 


- 1.28 12 


PBE 

Au(lll) OBS 
G06 
TS 


3.50 1.30 
3.48 1.30 

3.28 1.30 


2.31 1.29 18 
2.24 1.30 19 

2.23 1.30 18 


PBE 

Ag(lll) OBS 
G06 
TS 


3.64 1.30 
3.60 1.30 
2.75 1.33 
2.98 1.31 1 


2.27 1.32 23 
2.20 1.32 25 
2.14 1.32 25 
2.16 1.32 25 


PBE 

Cu(lll) OBS 
G06 
TS 


1.98 1.40 39 
1.97 1.40 38 
2.05 1.40 11 
2.05 1.40 13 


1.93 1.35 33 
1.91 1.35 33 
1.89 1.35 35 
1.89 1.35 34 



result unchanged by the dispersion corrections despite a 
small downward shift similar in magnitude to the changes 
of z observed at Au(lll). For the trans isomer, correc- 
tion effects range from none in the OBS scheme, to a dra- 
matic reduction of z with concomitant elongation of c^nn 
in the G06 and TS schemes. Nevertheless, the trans iso- 
mer at Ag(lll) again remains planar in all cases. Finally, 
at Cu(lll) the role of dispersion corrections is decidedly 
different from that at Au(lll) and Ag(lll): With GGA- 
PBE already yielding a strong azo-bridge to surface bond, 
which significantly elongates the NN bond and pins the 
bridge at a short vertical distance to the surface, the 
cis isomer z-shift induced by the dispersion schemes is 
smaller than at the other two substrates. Since this azo- 
bridge bond also determines the surface distance in the 
trans isomer and thereby puts the phenyl-rings well inside 
the range of surface Pauli repulsion, the intra-molecular 
distortion energy succumbs and as indicated by the large 
value for w, the phenyl-rings are bent out of the molec- 
ular plane in GGA-PBE. Again, the OBS scheme is too 
weak to influence this result. However, the G06 and TS 
schemes are not, and bend the phenyl-rings back towards 
the surface, yielding an w of some ten degrees, consider- 
ably closer to that of the gas-phase geometry (0°) than 
that of GGA-PBE (- 40°). In fact, this restoring force, 
self-consistent with the ring repulsion and molecular dis- 
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FIG. 3; (Color online) Schematic illustration of the geome- 
try changes induced by the G06 and TS correction schemes. 
Shown are sideviews with the GGA-PBE geometries shaded 
in the background as reference. The OBS scheme does not 
have a significant effect on the GGA-PBE geometries. 



tortion energy, is sufficient to somewhat offset the ef- 
fect of the azo-bridge— surface bond, which appears as a 
slightly increased surface distance z. 

As schematically summarized in Fig. [3] the OBS 
scheme correction is thus overall too weak to significantly 
modify molecular geometries. Given its shallow, long- 
ranged damping function, this result is readily under- 
stood: As apparent from Fig. [T]the less defined minimum 
in the correction potential turns the substrate dispersion 
interactions into a smooth background potential, with 
small gradient contributions. On the other hand, the 
GOG and TS schemes use similar, deeper damping func- 
tions with gradient corrections strong enough to make 
adsorbate geometries in principle dependent on the lo- 
cal substrate geometry. However, in the cis isomers the 
phenyl-rings sit largely outside the large-gradient range 
of this G06/TS damping function. The geometry correc- 
tions are therefore small and practically identical for the 
two schemes. For the planar trans azobenzene this situ- 
ation is different and the phenyl-rings do fall inside the 
large-gradient range. At Au(lll) and Ag(lll), where the 
weak azo-bridge— surface bond does not fix the molecular 
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TABLE III: Adsorption energies -Eads and work functions $ 
for the azobenzene cis (C) and trans (T) isomers at the three 
coinage metal surfaces, as well as the relative C-T energetic 
stability AE. Note that the G06 scheme does not feature 
parameters for Au. All numbers in eV. 







Eads 




$ 






T 


n 


C-T 


T 


n 




PBE 






0.57 






Gas-Phase 


OBS 






0.44 








G06 






0.47 








TS 






0.49 








PBE 


-0.12 


n 97 


0.42 


4.82 


o.oy 


Au(lll) 


OBS 


1 


u. y y 


0.50 


4.82 


o.oy 




G06 














TS 


-1.67 


-1.23 


0.93 


4.70 


3.89 




PBE 


-0.11 


-0.42 


0.26 


4.20 


3.68 


Ag(lll) 


OBS 


-1.05 


-1.18 


0.31 


4.18 


3.72 




G06 


-2.20 


-1.58 


1.09 


4.05 


3.73 




TS 


-1.76 


-1.41 


0.85 


4.05 


3.73 




PBE 


-0.27 


-1.08 


-0.24 


4.13 


3.87 


Cu(lll) 


OBS 


-1.97 


-2.25 


0.16 


4.12 


3.86 




G06 


-3.00 


-2.30 


1.18 


3.81 


3.87 




TS 


-2.81 


-2.46 


0.85 


3.81 


3.88 



height, the resulting attractive dispersion interaction cor- 
respondingly pulls down the entire molecule. At Cu(lll), 
the strong azo-bridge— surface bond prevents this, and 
the dispersive attraction only flattens the molecule by 
bringing the phcnyl-rings closer to the surface. 



C. Electronic structure and adsorbate energetics 



Despite these partly severe changes of the molecular 
adsorption geometry, it is interesting to note that central 
electronic structure quantities like the work function $ 
or projected frontier orbital positions are not much af- 
fected. Compared to corresponding GGA-PBE results 
detailed in our preceding study™ computed work func- 
tions in dispersion corrected geometries differ e.g. be- 
tween ~ 0.1 — 0.3 eV, while the entire qualitative ratio- 
nalization in terms of hybridization and charge rearrange- 
ment remains untouched. For completeness we therefore 
only present the quantitative work function values in Ta- 
ble mil but refer to the preceding publication for further 
details on the surface electronic structure^. In contrast, 
and not surprisingly, the adsorption energetics also sum- 
marized in Table [mi are substantially affected by the dis- 



persion correction schemes. The general pattern of these 
correction effects follows that observed for the adsorption 
geometries: At Au(lll), the pure GGA-PBE adsorp- 
tion energies are corrected by approximately 1 eV, with 
a comparable correction for both isomers with OBS, and 
the trans isomer correction some 30% higher in the TS 
scheme. Consequentially, the GGA-PBE relative isomer 
stability AE' is practically unaltered in the OBS scheme, 
but is doubled in the TS scheme. The Ag(lll) results 
are qualitatively as well as quantitatively comparable to 
those at Au(lll): Whereas OBS does not differentiate 
significantly between isomers, G06 and TS do, binding 
the trans isomer stronger and correspondingly increasing 
/S.E. Also at Cu(lll) the dispersion corrections quali- 
tatively resemble those at the previous two substrates, 
but are about 1 eV larger in magnitude. They are there- 
with in all three schemes so large that they restore the 
energetic order of the two isomers back to that in the 
gas-phase, i.e. the inversion with a more stable cis iso- 
mer obtained with GGA-PBE does not prevail. 

Apart from dramatically increasing the adsorption en- 
ergies, the overall prediction of the three semi-empirical 
schemes compared to the pure GGA-PBE numbers is 
thus a larger stabilization of the trans isomer upon ad- 
sorption. Within the GOG and TS scheme this dispersion 
corrected over-stabilization of the trans isomer is in fact 
so large, that the relative cis-trans energetic stability IS.E 
at all three substrates is larger than the corresponding 
gas-phase value. Despite this qualitative agreement, the 
corrected energetics provided by the three schemes shows 
quite pronounced quantitative differences. This holds in 
particular for the OBS scheme, which in the worst cases 
yields adsorption energies that differ by almost 1 eV from 
those of the G06 or TS scheme. The latter two schemes, 
on the other hand, are sometimes even more consistent 
in their results than could be expected from the variation 
of their corresponding scaled sq Ce^ij coefficients in Table 
m Considering that the variation to the OBS coefficients 
is of similar magnitude, we conclude that the decisive 
factor for the large difference in the obtained energetics 
from OBS vs. the G06/TS schemes is less the scatter 
in the dispersion coefficients, but the differences in the 
damping function. 



IV. DISCUSSION 

In light of the preceding first-principles GGA-PBE re- 
sults on the nature of the azobenzene— surface bond^'3. 
the intuitive effects one would expect from an additional 
attractive dispersion interaction on the adsorption ge- 
ometry and energetics are easily summarized: Adsorp- 
tion energies should become more exothermic, some- 
what more for the trans isomer with its phenyl-rings at 
closer distance to the surface. The additional interaction 
should pull the molecule further down, or, where this is 
prevented by strong covalent azo-bridge— surface bonds, 
should at least tend to flatten the molecular geometry 
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by driving the Pauli-repelled phenyl-rings to smaller z- 
heights. Per construction, this is exactly what the three 
semi- empirical schemes do when applied to the azoben- 
zene at coinage metal surface problem, and in this re- 
spect the results appear plausible. However, much more 
conspicuous than the qualitative trend is the sheer mag- 
nitude predicted. This is most pronounced for the ad- 
sorption energies, which in some cases are increased by 
more than 2eV. Despite the intended nature as a semi- 
empirical 'correction', the dispersion schemes thereby ac- 
tually provide a contribution to the final adsorption en- 
ergy that is up to an order of magnitude larger than the 
original first-principles result. 

Such results can not be uncritically accepted, in partic- 
ular when recalling the initially summarized fundamental 
doubts on the applicability of these correction schemes 
to bonding at metal surfaces at all. In the absence of 
detailed structural and energetic data for the azoben- 
zene problem either from experiment or high-level the- 
ory and as the phenyl-rings play the main part in the 
trends described, we exploit the similarity of the phenyl- 
ring moiety (— CgHs) to the benzene molecule (CgHg) to 
get at least some guidance on the order of magnitude 
for the vdW interactions. For the latter molecule TPD 
experiments indicate an adsorption energy of ^ 0.6 eV 
at all three substrates, albeit with rather large error 
barsJ^ii^'i^. Obtaining again only a minor lateral corru- 
gation of the potential energy surface, we focus here on 
the so-called 'HCP-B' benzene adsorption geometryi^, in 
which the molecule is parallel to the surface^, centered 
over the hep three-fold hollow site, and rotated such that 
three carbon atoms are maximally coordinated to the 
three-fold hollow metal atoms. 

For this geometry. Fig. 3] shows the computed binding 
energy curves for DFT LDA and GGA-PBE, as well as 
for the three semi-empirical correction schemes. In these 
calculations only the vertical z-height of the molecule was 
constrained, while all other degrees of freedom were fully 
relaxed. As in the azobenzene case, the semi-empirical 
schemes yield sizable corrections at all three substrates 
that increase the essentially zero GGA-PBE adsorption 
energies to roughly the spurious LDA value. Also, the 
trend predicted by the three dispersion schemes is very 
similar, with OBS adsorption energies consistently lower 
than the G06 and TS ones. In fact, compared to the data 
for the more similar trans azobenzene isomer in Table [III 
the correspondence is almost quantitative: For instance, 
at Au(lll) the adsorption energy correction with respect 
to GGA-PBE induced by the OBS and TS schemes is 
0.46 eV and 0.72 eV respectively for benzene, while it is 
0.93 eV and 1.55 eV for azobenzene with its two phenyl- 
rings. 



This corroborates benzene adsorption as a suitable ref- 
erence system and it is therefore meaningful to assess 
the performance of the three semi-empirical dispersion 
schemes with respect to the benchmark provided by the 




1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 



FIG. 4: (Color online) Binding energy curves of benzene, as a 
function of the center-of-mass surface distance z, at all three 
substrates. Shown is data for LDA and GGA-PBE, as well as 
the three semi-empirical dispersion correction schemes evalu- 
ated at ^-intervals of 0.25 A, to which the curves have been 
fitted as a guide to the eye. For each substrate, the exper- 
imentally measured TPD desorption energy is marked as a 
grey dashed horizontal line (see text). 



TPD experiments. For a more quantitative comparison 
Table IIVI compiles the corresponding data of the minima 
of the binding energy curves and additionally accounts 
for zero-point energy shifts. As obvious from these re- 
sults, all three schemes certainly provide an improvement 
compared to the GGA-PBE data, with the OBS scheme 
even falling close to the experimentally measured value at 
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TABLE IV: Results for benzene at the various substrates, in- 
cluding the experimentally measured adsorption energy , 
the computed adsorption energy without -Eads and with zero- 
point energy correction E^^^, as well as the optimum molec- 
ular height z. 









-Bads 


pZPE 
-L/ads 


z 




(eV) 




(eV) 


(eV) 


(A) 






PBE 


-0.08 


-0.03 


3.77 






LDA 


-0.73 


-0.68 


2.75 


Au(lll) 


-0.64i^ 


OBS 


-0.54 


-0.50 


3.48 






G06 












TS 


-0.80 


-0.77 


3.25 






PBE 


-0.06 


0.00 


3.75 






LDA 


-0.72 


-0.67 


2.66 


Ag(lll) 


-0.58" 


OBS 


-0.55 


-0.48 


3.38 






G06 


-0.97 


-0.92 


2.80 






TS 


-0.80 


-0.75 


3.01 






PBE 


-0.06 


-0.01 


3.48 






LDA 


-0.96 


-0.90 


2.18 


Cu(lll) 


-o.eaii 


OBS 


-0.82 


-0.75 


3.03 






G06 


-1.00 


-0.92 


2.63 






TS 


-1.05 


-0.98 


2.75 



all three surfaces. The other two schemes exhibit instead 
a tendency to overbind which gets most pronounced at 
Cu(lll). At this surface this in fact amounts to an over- 
correction by ~ 0.3 — 0.4 eV particularly by these two 
presumably more refined correction schemes. Within the 
simple analogy considering only the number of phenyl- 
rings, this would suggest that the corresponding trans 
azobenzene at Cu(lll) adsorption energies in Table |TT] 
could contain sizable errors of the order of ^ 0.6 — 0.9 eV. 

This trend of increased overbinding for benzene at 
Au(lll) to Cu(lll) goes hand in hand with a sys- 
tematically smaller height z of the molecule above the 
surface, cf. Table IIVI This brings the corresponding 
molecule— metal atom distances entering Eq. ([T]) closer 
and closer to the sum of their vdW radii and therewith 
into the distance range most affected by the damping 
function, cf. Fig. [T] In order to assess the implications 
also for the azobenzene problem, we note that at these 
intermediate distances this dependence on f{Rij , Rij) ex- 
tends not only to the actual functional form employed, 
but also sensitively to the cut-off radii^i^. Thus, while 
e.g. the metal-carbon correction energy at long range 
varies exclusively with the rather well motivated Cg^y 
coefficients (listed in Table |T| , it may easily in- or de- 
crease by a factor two with just a ± 10% variation in i?° 
when the metal-carbon distance is comparable to i?^^. 
We would therefore again attribute most of the observed 



scatter among the three semi-empirical schemes to the 
specifics of the employed damping function. 

In this respect it is discomforting to realize that this in- 
ffuential part of the semi-empirical correction expression, 
Eq. ID), is the one least motivated from first-principles. 
Introduced to merely prevent double counting of those 
dispersion interaction contributions that are already con- 
tained in the employed local or semi-local DFT xc func- 
tional at short and intermediate range, there is no guid- 
ance on the exact form, nor even existence, of an ana- 
lytical function that would achieve this. The difference 
in our OBS versus the GOG and TS results shows that a 
certain steepness of the damping function is required to 
achieve a significant dispersion-induced change of the ad- 
sorption geometry. On the other hand, we verified that 
the dependence on the onset of the damping function, 
i.e. the cutoff-radii, is intuitive. Larger cutoff-radii re- 
duce the dispersion correction, which means for the sys- 
tems studied here that the adsorption energy is reduced 
and the adsorption geometry comes closer to the origi- 
nal GGA-PBE one. However, without knowing the exact 
structure and binding energy, it is impossible to conclude 
in which direction the damping function employed by the 
different schemes would need to be modified, even if other 
uncertainties of the approach like the neglect of metallic 
screening could be excluded as the source of the error. 

Only corresponding detailed data from experiment or 
high-level theory can therefore provide the final answer 
as to the performance of the semi-empirical correction 
schemes for the azobenzene at coinage metal surface 
problem. In the absence of such data we cautiously 
conclude from the comparison to the similar benzene 
adsorption problem that the three schemes employed 
in this work seem to provide an account of dispersive 
vdW interactions that is in the right ballpark, and could 
even be semi-quantitative. The OBS scheme is likely 
mostly limited by its shallow damping function that 
yields too small gradients to significantly modify the ad- 
sorption geometries. GOG and TS on the other hand are 
likely to overbind, with a maximum estimated error of 
~ O.G — 0.9 eV in the adsorption energies. 



V. SUMMARY 

The problem of azobenzene adsorption at coinage met- 
als combines a large molecule, heavy transition metal 
surface chemistry and a multi-facetted bonding mecha- 
nism, in which dispersive vdW interactions play a crucial 
role. This combination forms a tremendous challenge for 
contemporary first-principles modeling, with DFT with 
semi-local xc-functionals still setting the standard. In 
order to assess the effect of the insufficient description 
of dispersive interactions at this level we employed three 
different semi-empirical correction schemes that account 
for such interactions at least on the level of the lead- 
ing long-range term. The low computational cost 
of these schemes makes them an attractive solution for 
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this problem, and the overall trend obtained by applying 
them to the azobenzene problem at Cu(lll), Ag(lll) 
and Au(lll) is consistent with the anticipated effects of 
an additional attractive interaction: Compared to pure 
GGA-PBE adsorption energies become more exothermic, 
somewhat more for the trans isomer with its phenyl-rings 
at closer distance to the surface and thereby reinforcing 
the higher gas-phase trans isomer stability. The addi- 
tional interaction pulls the molecule further towards the 
surface, or where this is prevented by strong covalent 
azo-bridge— surface bonds, at least flattens the molecu- 
lar geometry by driving the Pauli-repelled phenyl-rings 
to smaller heights. 

Much more problematic than this overall trend is the 
sheer magnitude of the effects predicted. Constructed 
as a semi-empirical 'correction', the dispersion schemes 
thereby actually provide a contribution to the final ad- 
sorption energy that is in some cases up to an order 
of magnitude larger than the original first-principles re- 
sult. In addition the three conceptually similarly rooted 
schemes exhibit a discomforting scatter in their results 
that amounts up to 1 eV in the adsorption energies and 
up to 1 A in central geometry parameters. An analysis of 



the congeneric adsorption of benzene at the three coinage 
metal surfaces suggests that the three schemes employed 
in this work seem to provide an account of dispersive 
vdW interactions in the azobenzene adsorption that is 
in the right ballpark, with the more refined G06 and TS 
schemes likely to overbind. The ultimate answer as to the 
performance can only be provided by comparison to de- 
tailed structural and energetic data from high-level the- 
ory or experiment, which could then critically stimulate 
further development and improvement of these schemes. 
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